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ABSTRACT: We consider here renormalizable theories without relevant 
couplings and present an I.R. consistent technique to study corrections to 
short distance behavior (Wilson O.P.E. coefficients) due to a relevant per- 
turbation. Our method is the result of a complete reformulation of recent 
works on the field, and is characterized by a more orthodox treatment of U.V. 
divergences that allows for simpler formulae and consequently an explicit all 
order (regularization invariant) I.R. finitess proof. Underlying hypotheses 
are discussed in detail and found to be satisfied in conformal theories that 
constitute a natural field of application of this approach. 
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1 Introduction 



We present here a consistent technique to study the corrections to the 
short distance behavior of a renormahzable euchdean quantum field theory 
without relevant couplings (i.e. with positive mass dimension), when a rele- 
vant perturbation is introduced. 

It is well known that expansions in relevant parameters around field the- 
ories without a mass scale are plagued by I.R. divergences. Considerable 
effort on how to deal with these singularities has been done, for generic 
super-renormalizable quantum field theories [0 as well as for some specific 
cases [|], but the present situation is such that any general and (especially) 
practical method that could face these I.R. divergences is always welcome. 

The possibility of doing an I.R. finite expansion in a relevant parameter 
is particularly appealing in the case of (two dimensional) conformal field 
theories. After the first pioneering work [Q, understanding on this subject 
has been greatly increased and many exact results have been found (see for 
instance |Q). In addition it soon became clear that those models can be 
taken as a starting point to get information on non conformally invariant 
theories 0- From a point of view of statistical mechanics, 0, conformal 
field theories are fixed points of the renormalization group and the line of 
research we are considering can be rephrased as "how to go outside the fixed 
point" . We think that the technique introduced here finds a natural field of 
application within this framework because it gives a general way to describe 
the short distance behavior of a theory outside the fixed point, useful to test 
eventual exactly solved deformations but essential to give new light in all the 
other cases. 

The method described in this paper is a reformulation of the independent 
works of Al.B. Zamolodchikov and of Mikhak and Zarkesh (that 
applied techniques developed by Sonoda 0). 

The first idea underlying the approach is present already in earlier papers 
of Wilson [jll| and Wegner |]12| (see also the more recent [|13[)- to deal with 



I.R. singularities is essential to redefine composite fields in such a way to 
confine the terms non analytic in the generalized masses (potentially present 
in all the interesting correlators) inside the V.E.V. of the composite fields, 
enforcing in this way the regularity of Wilson coefficients. 

The second step is to give a formal expression of the derivatives of Wil- 
son coefficients with respect to the generalized masses in terms of (generally 
unknown) correlators of the deformed theory: this has been realized by per- 
turbative expansions in 0], and by using the "Variational Formulae" in ^ 
(a kind of Action Principle, see Section El). 
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The third step is to regularize the theory with an I.R. cutoff and to take 
the hmit of zero generahzed masses: clearly at this stage one has to face the 
potential I.R. divergences. The other main idea of the approach, [||, ^ (and 



10 1 for an explicit proof up to second order), is that the convergence prop- 
erties of O.P.E. should guarantee the cancellation of I.R. divergences. After 
keeping into account this fact the derivatives of the Wilson coefficients can 
be computed by evaluation of I.R. cutoff integrals involving exactly known 
correlators of the massless field theory. 

Our present contribution is an improvement and an extension of the pre- 
viously mentioned works, in order to give a simplified and U.V. consistent 
method with well defined applicability conditions that is I.R. finite at all 
orders. 

First in Section ^ we state clearly which are the underlying, minimal 
hypotheses to be satisfied for the method to work. 

Secondly, by using the orthodox renormalized Action Principle (that has 
a simpler form with respect to 0), an explicit inductive proof of I.R. finitess 
and a detailed discussion of the implications of cancellation of I.R. diver- 
gences is given in Section ^. It is also shown that final expressions do not 
depend on the concrete form of the I.R. regularization (whose choice becomes 
a matter of convenience). 

In Section ^ we treat the correct definition of renormalized composite 
operators to be chosen in such a way that Action Principle holds. We shall 
see that from this requirement arises the problem of existence of numerical 
constants (essentially present in the operators' V.E.V.) whose expression in 
terms of deformed theory correlators is well defined but in general uncom- 
putable without additional information. 

In Section we compute the mass corrections to the short distance be- 
havior of the propagator of two dimensional free euclidean fermions, while in 



Section |5.2| we study as a further example of our approach the mass correc- 



tions to the spin spin correlator in the Ising model. 



2 The general method and underlying hypotheses 

We present in this Section the ideas of the approach (see also ||, |, p!0| ) 
and state a clear list of the underlying hypotheses (which are essentially 
three) . 

The goal is the reconstruction of the short distance behavior of corre- 
lation functions of composite operators of a D dimensional euclidean field 
theory without dimensional couplings, perturbed by one or more operator 
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Oi (of canonical dimension Xi) with dimensionful couplings m* ("generalized 
masses" of dimension yi = D — Xi): 



AS = - jdxY^m'sOU^), (2.1) 

i 

(notice the minus sign above, taken for convenience). We will refer to 
this model as "deformed theory" while the unperturbed one will be called 
"massless theory". The perturbation operator Oi is supposed here relevant 
(0 < Xj < D) to restrict the number of renormalization conditions needed to 
define its multiple insertions, see Section |[ 

We will assume that all correlators of the deformed and of the massless 
theory to be well defined. 

The massless theory is supposed to be at least perturbatively U.V. renor- 
malizable with respect to some dimensionless coupling A; furthemore we will 
require that only logarithmic corrections to tree level scale invariance can 
arise. 

The eventual perturbative (with respect to A) knowledge of correlators at 
m* = gives a perturbative estimation of each term of the mass expansions 
of the correlators of the deformed theory; nevertheless in the particular case 
of (two dimensional) conformal field theories (bare) correlators are exactly 
known and thus coefficients of mass expansion can be exactly computed. In 
the following we will forget the eventual dependence in A with the convention 
that all formulae should be intended to hold order by order in renormalized 
perturbative expansion with respect to A in the general case. Notice that the 
absence of super-renormalizable couplings guarantees that I.R. problems do 
not arise in massless theory if a renormalization scheme without additional 
unphysical masses is used (see e.g. discussion in the second and third papers 
of 1T§). 



The short distance behavior of the deformed theory is described by the 
Operator Product Expansion ||11|| : 



< [<I)^,(ri) (r„)]X Q...^Jri-r„,---r„_i-r„;m) < [$c(r„)]X >„ 

(2.2) 

where are a complete set of composite operators (that reduce in the limit 
m* — > to the operators of the massless theory), X is a multilocal product of 
other fields localized "far away" from the r^, the suffix m refers to deformed 
theory, and the convergence of the series will be discussed below. In the 
following we will denote the renormalized multiple insertions of operators, 
see Section H, enclosing them with [•■■]. 

We will present an I.R. finite technique to compute multiple derivatives 
of Wilson coefficients C^^...„^ with respect to (renormalized) m* at the point 
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m* = in terms of integrated correlators of the massless theory. We will 
assume the following: 

Hypothesis 1 (Regularity): an U.V. renormalization scheme for cor- 
relators of the deformed theory is assumed such that counterterms are polyno- 
mial in renormalized generalized masses (Minimal Mass Dependence, MMD). 
Furthermore, the Wilson coefficients are supposed to be regular (i.e. C°°) in 
generalized masses at m' = 0. 

Minimal Mass Dependence guarantees the smoothness of the m* — > 
limit from the point of view of U.V. renormalization, property which will 
be essential in the following. Any subtraction scheme that subtracts only 
the infinite part in the underlying regularization (Minimal Subtraction) is 
expected to satisfy MMD because no additional mass dependence is intro- 
duced in counterterms (apart from the trivial powerlike factors) that gives 
the correct physical dimension. 

We remark that Wilson coefficients and Minimal Mass Dependence are 
not independent: at least in the framework of perturbative renormalization of 
quantum field theories MMD is equivalent to regularity of Wilson coefficients. 
This can be understood as follows: first we remind that Minimal Subtraction 
in Dimensional Renormalization (MS) satisfies MMD, and regularity, 
2D[. Secondly any new scheme satisfying MMD will differ from MS by means 



of finite counterterms, polynomial in the masses. In particular the change of 
scheme for renormalized operators will amount to the transformation 

[*a]' = (2.3) 

with Nj!^ = (5^ + and Y a matrix of order 0{h) that is polynomial in the 
masses. Wilson coefficients will change accordingly as 

C'l...a„ = Kl ■ ■ ■ <"iV-^^Ct..,„ (2.4) 

and the new ones will be clearly regular in m = as the ones in MS scheme. 
Conversely in a scheme with non regular counterterms the matrix A*"^ above 
will have some non regular entry and so for Wilson coefficients. 

We want to emphasize furthermore that the stronger (nonperturbative) 
assumption of analyticity of Wilson coefficients in generalized masses would 
ensure the convergence of the Taylor series we are building up here, which 
is only asymptotic in the general case. 

The second assumption necessary to compute derivatives with respect to 
generalized masses is: 

Hypothesis 2 (Action Principle): it is assumed that for each renor- 
malized generalized mass nf a conjugate operator Oi exists such that the 
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derivative with respect to rrt is realized as: 

d^. < [X] >„= jdx < [: 0,{x) : X] : := O,- < [O,] >„ (2.5) 

for each multilocal operator X. 



The Action Principle [|IJ] is well known to hold in perturbative renormal- 



ization, in the classical B.P.H.Z. framework [15] as well as in other schemes 



satisfying Hypothesis 1 such as Dimensional Renormalization ||T6[ and An- 



alytic Renormalization |]T^ (see also |TH] for Differential Renormalization). 



We will come back to this argument in Section ^. 

The simplicity of Eq.( p^.5D compared to "Variational Formulae" |^ is the 
key of our approach that will allow us to deal with expressions for derivatives 
of generic order and to give the inductive IR finitess proof. 

Let us finally have a glance at the method in a simple case, to introduce 
the final hypothesis. Suppose we want to compute the first derivative with 
respect to m* (in m* = 0) of the two point Wilson coefficient (with a, 6, 1 
referring to renormalized operators 1). Let us consider the quantity 

(remainder of O.P.E.) 

Xc<N 

A(^)(XH,m) [($.(r)$,(0) - £ C:,$,(0))X^] >^ (2.6) 

c 

in which Xc is the canonical dimension of operator $c and Xr is the trivial 
operator {Xr = 1) or is a multilocal operator with support in 

Er = {x/\x\ > R}. (2.7) 

such that R > \r\. (In the following we also denote by Ir the complement 
set oiER, i.e. Ir = R''-Er). 

We then take derivatives of both sides with respect to m*: by using Action 
Principle we can write 

„ Xc<N 

dA^''\XR,m) = dx < [: 0,{x) : ($,$, - ^ C'^Xr] >^ 

•' c 

Xc<N 

- E ^^C:^b < ['^cXr] >™ . (2.8) 

c 

It is evident that the in the limit of zero generalized masses the integral 
in the right hand side might be I.R. divergent. To deal with this problem we 
must proceed as follows, 0. First we split the integral in two pieces, 

9,AW(X^,m) = / dx<[:Or.{<l>a<^b- C:MXR]>m 
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+ / dx<[.Or. ($a$6 - E C:,^c)Xr] >^ 

Xc<N 

- £ d,C:, < [^,Xr] (2.9) 

c 

in which Ri > R. Then we observe that the second integral can be rewritten 

as 

/ dx<[:Or. ($a$6 - E C:,<I>c)Xr] >m= AW(X;j,m) (2.10) 

i.e. in the same form of Eq. (|2.6| ) with X'^ = Ie '■ '■ satisfying the 
same hypothesis on the support as Xj^. 

If we could assume the (weak) convergence of O.P.E. series when inserted 
in correlation functions containing "far" operators Xji for the deformed the- 
ory, we would say that the limit oo of A^'^\Xji,m) of (|2.6|) is zero 
together with all its mass derivatives and in particular the N —>■ oo limit of 
Eq.( p.lOD would be zero for arbitrary i?i > i? > |r|: the subsequent limit 
m — > of Eq. (|2.9|) would be safe and could be exchanged with the first (I.R. 
cut off) integral with the desired result 

„ Xc<N 

= lim { / dx<[:Or. ($a$6 - E C:MXn\ >„=o 

Xc<N 

- E diC:, < >„=o} (2.11) 

c 

(the suffix <>^=o will be omitted in the following). 



We think that, in spite of some argument in favor |£T|, convergence of 
O.P.E. is far from being reasonably proved in a quantum field theory with 
(generalized) mass; in particular the existence of terms of the form e~^^^"^'^^^ 
could spoil the convergence of the expansion for small r (see the ffist paper 
in i). 

Nevertheless for our purposes it suffices the following weaker hypothesis 
that we state here for general Wilson coefficients. 

Hypothesis 3 (O.P.E. asymptotic weak convergence): it is as- 
sumed that the remainder of O.P.E. 

Xc<N 

AW(X^,m)=<[($.,(rO (r„)- £ Q...a„<fc)X«] >rn (2.12) 
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(X[i being the unity or an arbitrary multilocal operator with support outside 
R> |ri|, ■ ■ ■ \rn\) satisfies 

lim lim lim di, ■ ■ ■ ^^,A^^\XR,m) = (2.13) 

iV— >cxD m^O R^oo 

for each k. Equivalently [| we can write that 

Xc<N 

lim lim < [($„,(ri) ■ ■ ■ $,Jr„) - £ Q...a„<f c)^^] >~ (2.14) 



in the sense of asymptotic series (in nf). 

First of all we notice that assumption Eq. (|2.14| ) for = requires the 
convergence of the O.P.E. in the underlying m* = theory when inserted in 
correlators with "far" operators. This property is well known to hold in con- 
formal field theories, WM- Also from axiomatic PSI as well from perturbative 



24| considerations we can say that in general quantum field theories O.P.E. is 
an asymptotic expansion in powers of r. But in the limit R oo, by cluster- 
ing, the contribution of in the correlator factorizes and the only physical 
scale in the O.P.E. becomes the generalized mass (the renormalization point 
IX cancels between V.E.V. and Wilson coefficients by Renormalisation Group 
invariance of the sum) and a remainder of order 0{r^) is actually of the form 
0{{mvry) by dimensional analysis, and thus the O.P.E series with respect 
to powers of the generalized masses reasonably becomes asymptotic. 

Let us come back to our previous example. In the limit of Eq.( p.l3 ) it is 



clear that the contribution from Eq. (p.lO|) vanishes. Observing that also diA 
vanishes, and choosing = 1 we obtain the formula 



Xc<N Xc<N 

= lim / <[:Or. (<f a<f E c:,m,)] > - £ d.c:,{o) < > 

(2.15) 

in which the limit over m and R were exchanged without problems for the sur- 
viving I.R. regulated quantities. (Notice also that (0) in Wilson coefficients 
refers to m = 0; we will omit this specification for simplicity of notations in 
next Sections, being clear from the context if we are dealing with zero mass 
Wilson coefficients or not.) 

To have a useful formula we must observe that in general, by dimensional 
considerations, 

lim / rfxi ■ ■ ■ / dxk< [: Oil : ■ ■ ■ : O., : $^(0)] >= (2.16) 



^ In the rest of the paper we will exchange derivatives with respect to generalized 
masses with the limit over I.R. cutoff when the limit over m follows, because at to 7^ the 
limit over R is smooth. 
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if Xc — Z]j=i Vj > {Vi being mass dimensions of m*) due to the absence of a 
physical scale in the m = theory (we are using the fact that renormalization 
point /i does not give powerlike corrections). 

By Eq. ( p.l6D the series in c index is truncated and we can write 



^c<yi 



d,Cl,{0) = hm / dx<[:Or. ($a$6 - E C:,mc)] > . (2.17) 

The goal of expressing the (first) derivative of a Wilson coefficient in term 
of well defined and known quantities of massless theory is then reached! In 
the next Section we will deal with higher order derivatives. Notice that 
all formulae that we will derive do not depend on the U.V. renormalization 
scheme, provided that Hypotheses 1-3 are satisfied. The choice of the scheme 
is a matter of taste and computational convenience. 



3 All order formulae and I.R. finitess proof 

In this Section we will give I.R. safe expressions for the n^^ derivative 
of the Wilson coefficients with respect to generalized masses by using the 
assumptions of previous Section. Then some consequences of these identities 
will be discussed. 

We will prove by induction on n the general form of ra*^ derivative of (two 
operators) Wilson coefficient at zero masses and its I.R. finitess. Generaliza- 
tion to higher operator Wilson coefficients is straightforward. 

We want to prove that at order n the following relation holds for the 
deformed theory: 

lim lim 9,,---9,„A(^)(XR,m) ~ 



x,,<N 



lim { /dxi ■ ■ ■ / dxn < [: : ■ ■ ■ : O., : ($.,$., - £ C7i„^$,)X«] >, 

- E ^nCaiaa / ■ ■ ■ / < [: : ■ ■ ■ I Oi^ : ^^Xr] >^ +p. 

J I2 J In 

- E 9iA2Cl,a2 dx3--- dxn<[:0,^:---: : ^^Xr] >„ +p. 

J I-t J In 



xi,<N 



E ^.r■■^^nCla2<[^bXn]>m} (3.1) 

in which = Ir. for shortness and Rn-i ■ ■ ■ -Ri > -R > |r|. The series are 
to be intended asymptotic as explained in the previous Section and with p. 
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we mean all terms obtained interchanging the suffix (1, ■ ■ ■ , /c — 1) of Wilson 
coefficients' derivatives with {k, ■ ■ ■ ,n) inside integrals without distinguishing 
the ordering. We observe that the left hand side of the previous expression 
is zero asymptotically due to Hypothesis 3, Eq. (|2.14| ), and that due to the 



presence of I.R. cutoff the massless limit is I.R. safe as we shall use below. 
It is clear by Hypothesis 3 that at n = Eq.(|3.1|) holds. Now let us 



assume that Eq.(3J-) holds at order n and show that it holds at order n + 1. 
Let us take derivative with respect to m*"+i of both sides of Eq. ( |3.1[ ) (and 
commute derivative with respect to limit over R, see footnote |l|). 

The action of i9j„_|_i on the correlators gives terms of the form 



lim 



jdxn+i < [. Ci„+i(a;„+i) : {n}] >m (3.2) 



i.e. the same terms of order n (contracted with Wilson coefficients as in ( |3.1| )), 
indicated above with {n}, but with the additional insertion of : Oi„+i :, and 
terms in which derivative acts on Wilson coefficients present at order n. But 
by Eq.( ^) (with Xr Jj^^^^dxn+i : Ci„+i : Xr) we obtain 



N. 



lim dxn+i < [■■ Oi„^, : {n}] lim / dxn+i < [■ Oi^^, : {n}] >^ 

(3.3) 

if we restrict to -R„-i, ■ ■ ■ , Ri > R ^ oo. 

It is easy to check that the sum of the two contributions simply reproduces 
the n + 1 expression according to (|3.1|) so that induction works Q.E.D. 

If we choose Xr = 1 in Eq. (|3.1| ), take the —>■ limit and remind 
the dimensional selection rule Eq. ( |2.16| ) we obtain finally that, in the limit 



Rn, Rn-i, ■ ■ ■ , Ri > R ^ oo, the series in index b are actually finite sums 
and that 

f). . . . f). — 

lim { / dxi ■ ■ ■ / dxn < [: : ■ ■ ■ : O., : ($.,$., - E Cl^a.^,)\ > 
-E^nC^ia. fdx2--- f dxr.<[:0,„:---: O,, : $,] > +p. 

^ •'^2 ''In 

- E 5n ■ ■ ■ d.^_,Cl,, [ dxn < [: : M > +p.} (3.4) 

b 

where < [: C'i,. : ■ ■ ■ : O^^ : $b] > is restricted to such Xb < Vi^ hl/i„. 
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More generally if we keep in Eg . (13.11) a nontrivial | and take the 
m* — i> limit, we can observe that 

I ■■■ I <[.0.„:---:0,,: $,(0)X^] >= 0{R') (3.5) 

•'Ik •'In 

in which we defined the I.R. degree of divergence of the correlator as 

n 

^ = E V'h -^b- xxji (3.6) 

j=k 

by simple dimensional considerations. It follows that all correlators with 
6 < give vanishing contribution in the R ^ oo limit and the series in the 
index b are effectively truncated. 

We obtain in all generality: 
= hm { / t/xi ■ ■ ■ / < [: : • • • : O,, : ($,,$,, - C'a,a2^b)Xn\ > 

- E dnCl^2 [dx,--- [ dx^<[:0,^:---: O,, : ^Xr] > +p. 

^ •! h •'In 

- E 5n ■ ■ ■ d.,^-A,2 I dx^ < [: : K.X^] > +p. 

b 

-j:dn---d^„C!:^,,<[<^bXn]>} (3.7) 

b 

where J^l < [• ^in ■ ■ ■ ■ ■ Oi^, : > is restricted to $b such Xf, < 

yik'-' + Vin -xxn- 

Eq. (|3.7| ) gives an infinite number of linear algebraic relations (one for each 
choice of Xr) among derivatives of Wilson coefficients up to order n. Of 
course consistency of the method (guaranteed by our assumptions) enforces 
that only a finite number of relations are independent and the system is not 
overconstrained. 

Moreover it is important to observe that in each equation such as ( p.7|) 
the I.R. divergences of the same form cancel independently. In particular this 

^ As an example, in conformal field theories the particular choice Xr = <i>c(-z, z)z^z^ 
(with z, z complex coordinates of modulus R and A, A the relative canonical dimensions 
of operator <i>c in the underlying conformal field theory) will select di^ - ■ ■ Si^C^^^j^ in the 
last term of Eq.( ^.7| ), by orthonormality property of conformal operators. 
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happens for terms of the same dimensionahty (I.R. degree of divergence): if 
for fixed ai , 02 we define the finite set 



k 

^ {d,, ■ ■ ■ d,,Cl^Ja = Y: + xt} (3i 



it is clear that for every choice of ii,---,in and Xr in Eq. (|3T7| ) we have 
linear relations only between elements of the same sets So-, because 6 = 



I]j=i Ui-j ~ ^Xr — o" for each term of (|3.7|) . 



At this point to compute a derivative of a Wilson coefficient belonging to 
a set 5*0- one has only to choose a convenient number of subrelations of ( p.7| ) 
(varying Xr or n) and solve a linear system. Notice that due to possible 
different form of I.R. divergences with the same dimensionality the number 
of relations obtained by I.R. cancellation might become even bigger. This 
can happen also when some quantum numbers are conserved. 

We want to show now that, while the particular I.R. spatial cutoff we 
used was essential to derive Eq.( p.l[ ) by induction, our final results, Eq.( ^^ 
and Eq.( p.7|) hold for a wide class of I.R. regularizations. The point is that 
in Eq. (|3.1| ) all the terms have the form (keeping only essential features) 

lim lim dxi ■ ■ ■ dxkOR^{xi) ■ ■ ■ 6R^^{xk) < ■ ■ ■ Xr > (3.9) 

m^O R.i>R^ Qo J J 

where 9r{x) = 9{\x\ — R) (with 9 the usual step function). Let us now 
consider a generic I.R. regulator function Qr/{x) {R' being a new length 
parameter) such that \imRi^oo^R'{x) = 1. Inserting 1 in the integrals and 
exchanging integration with these limits we have 

lim lim lim dxi--- / dxk9R^ (xi) ■■■9r^^ {xk)&R' (xi) ■ ■ ■ &r' (xk) < ■■■Xr> 

(3.10) 

It is clear that at m 7^ (the limit on m being external) we can exchange 
limits over R, Ri with those over R' and subsequently the limit over R with 
integrations, obtaining 

lim lim dxi ■■■ dxkQm (xi) ■ ■ • 0r' (x^) < ■ ■ ■ > . (3.11) 

Proceeding exactly as done previously (exchanging limit over masses and 
over R', using selection rules at m* = 0) we obtain for a generic I.R. cutoff 
function Qr: 



= lim { dxi--- dxn&R^ixi) ■ ■ ■ 6ij„(x„) x 

Ri^co J J 
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b 

h 

< [: C>,„ : ■ ■ ■ : : ^bXj\ > +p. + • • • 

b '' 

- E9,,---9,„Ci,^<[<l>,Xoo]>} (3.12) 

b 

(same meaning of the * as before). It is clear that a clever choice of I.R. 
cutoff can considerably simplify computations. For instance it is easy to 
prove that when X^o = 1 and the cutoff function 0r is rotation invariant, 
only scalar operators are involved in the expressions. This property will be 
used in Section ^ 

We close the Section noting that, while Eqs.(p.4D, (p.7|), were derived by 
using only of the condition Rn-i, ■ ■ ■ , Ri > R ^ oo, with the result of 
having symmetric formulae, for practical calculations it is better ^ to force 
the I.R. cutoff to go to infinity in a hierarchical way, say, Rn » Rn-i, ■ ■ ■ » 
Ri > R oo so that we can perform first the integral with bigger cutoff Rn, 
neglecting subleading terms of order 0{ ^^~^ ) and so on. This could not be 
the most convenient way to compute integrals in other I.R. regularizations: 
for instance in the case Gq^(x) = e**^*^ it is simpler to take all Qi = Q ^ 0. 



4 U.V. definition of composite operators 

In this Section we schematically review some known facts on the treat- 
ment of renormalized multiple insertions of composite operators with partic- 
ular reference to the operator Om conjugate to the generalized mass m by 
the Action Principle Eq.( p.5|) (we will deal here only with one generalized 
mass; extension to general case should be straightforward for the reader). 

The foundations of renormalization are locality and power counting (||25||, 
see p6| for a pedagogical introduction). Each composite operator has 
a related power counting dimension equal to the canonical one Xa, (we will 
not consider oversubtracted operators here). Each insertion of $a in a trun- 
cated connected one particle irreducible graph 7 changes the U.V. superficial 
degree of divergence by the amount 600 = —ya = —D + Xa, generating 
new divergences. Divergences related to one operator insertion are treated 
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allowing for a mixing of each bare operator with those of lower power count- 
ing dimension to get a finite renormalized one, Renormalized multiple 
insertions of composite operators (denoted here by [■■■]) contain new diver- 
gences that should be cured by counterterms proportional to delta function 
and its derivatives, the so called contact terms: for instance 

[$a(x)$fe(0)] = [<l>a(x)][$fc(0)] + contact terms. (4.1) 



For detailed treatment of arbitrary multiple insertions see and references 
therein. A practical way to deal with multiple insertions is obtained introduc- 
ing sources lj"'{x) of dimension ya for composite operators [^a], and allowing 
for counterterms proportional to monomial in the sources, in the composite 
fields and their derivatives of total dimension less then or equal D (see e.g. 



28|, m). 

For a composite operator to be well defined a normalization condition 
should be given in correspondence of each infinite subtraction. A renormal- 
ization scheme is a consistent and physical way to subtract divergences and 
to fix finite counterterms with explicit or implicit normalization conditions. 

Equivalently, a renormalization scheme can be seen as a consistent tech- 
nique to build up well defined tempered distributions (satisfying physical 
requirements, [^) from "bare expressions" that do not define a distribution 
and to fix finite counterterms arising from the ambiguities of this extension. 
This point of view is more suitable for conformal field theory where "bare" 
exact expressions are known for correlators of composite operators at differ- 
ent positions. 

For instance the bare correlator with two insertions of the operator £ = 
— : ipil) : in the free massless D = 2 euclidean Dirac Fermions model (that is 
a conformal field theory) is 

< S{r)S{0) >B= x^^^ (4.2) 

but 1/r^ is not a distribution. This can be understood because it is not 
locally integrable with respect to the measure cPr. A distribution T(r) that 
extends the bare correlator can be defined (in polar coordinates) as the so- 
lution of the distributional equation 

\r\T(r) = (4.3) 

{1/r is a distribution) whose general solution will have the form 

T{\r\) = Tp{\r\)+TH{\r\) (4.4) 
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where Tp is a particular solution of the complete equation such as -^P gTri^ff 
(Hadamard finite part, see [Appendix lj| ) and Th = Co6{\r\) is a solution of 



the homogeneous equation. The renormalized correlators will be equal to T 
and the dimensionless finite contact term Co will be fixed after imposing a 
normalization condition. The alternative definition 

|r|"T(r) = — |r|"-i (4.5) 

for some n > 1 would give a Th = Y.k=o Ck5^''\\r\) and would change dras- 
tically the short distance behavior of the correlator: in this sense is not 
physically relevant. 

From our discussion it appears clearly that only a finite number of in- 
sertions of a relevant operator give U.V. divergences (and contact terms to 
be eventually fixed): this is the main practical motivation to restrict to rele- 
vant perturbations, while all formulae presented in Section ^ can be extended 
formally to any marginal operator (at any finite order). 

After this brief digression let us come to the problem of defining renor- 
malized operators (i.e. of fixing schemes) satisfying the Action Principle. 
The point is that, while being clearly true at bare level (think about func- 
tional integration), the Action Principle might be not satisfied in a particular 
scheme, being possibly violated by renormalization. In this case the starting 
point scheme must be corrected (i.e. changed) by finite counterterms to en- 
force AP. It follows in our case that all the U.V. divergent insertions of the 
renormalized composite operator Om to be cured in a chosen scheme might 
require additional contact term to be fixed in such a way that Eq.(p.5|) holds. 
These contact terms result to be finite numbers (apart the trivial power of 
m) whose expression in terms of integrals of U.V. renormalized correlators 
(in the original scheme) is known and the only free parameter characterizing 
the composite operators is the usual renormalization point fi. 

Of course it is possible that in some renormalization scheme the Action 
Principle holds automatically without finite renormalization; this is the case 



of Dimensional Renormalization, where no additional contact terms are 
needed and life is easier. 

In our approach we are interested only in the zero mass limit, thus all 
dimensionful contact terms vanish and are not involved in our expressions. 
Only dimensionless contact terms could survive in zero mass correlators: 
they have a definite expression in terms of massive integrals but cannot be 
computed in terms of the massless quantities only. This situation arises when 
one tries to impose the condition 

9^ < >rn= JdXi--- dXk < [OrniXi) ■ ■ ■ 0„,{Xk)^a] >m (4.6) 
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in the case when Xa = ki/m (in general when the I.R. degree of divergence 5 
defined in Eq. (|3.6|) is zero). Equation ([4 .61) presents non trivial U.V. renor- 
malization and also is I.R. singular at m = (logm terms are clearly present 
in the operator V.E.V.): it follows that a numerical constant that parame- 
terizes the I.R. behavior of the operator in the deformed theory is left not 
fixed (but expressed in terms of deformed theory correlators). 

These unknown constants will parameterize our results. In the case of 
one generalized mass we can always confine this ambiguity to the composite 
operators V.E.V., as will be shown in the next Section with two examples. 
See also the related paper fS^ 



Before closing this Section we comment here that "Variational Formulae" 
of Ref. were invented with the goal of building up a geometrization 



of the theory space of the renormalized theories, |31|, with a connection 
expressed in terms of a set of finite arbitrary " contact terms" . We feel that 
the idea of geometrization (interesting by itself) is somewhat incompatible 
with straightforward calculations and in our pragmatic approach we will 
try to minimize the number of contact terms. In particular a locally flat 
coordinate system (no connection at all) is obtained when schemes with a 
built in Action Principle can be used. 



5 Two examples 

5.1 A toy model example: short distance behavior of 
two dimensional free fermion propagator 

In this subsection we compute the mass corrections to the short distance 
behavior of the two dimensional euclidean free Dirac fermions propagator 
up to order by using the knowledge of the massless theory. We will use 
as renormalization scheme the Minimal Subtraction (MS) in Dimensional 
Regolarization. 

To fix notations (we use that of [^) we write here the exact expression 

^ In our mind "Variational Formulae" are incomplete because do not give a general 
definite prescription to build up distributions from bare correlators (i.e. do not define a 
scheme). In particular it seems to us that there is no prescription to subtract singularities 
arising when three or more operators collapse to a point, as happens at high derivative 
orders. 
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for the propagator we want to recover: 

(5.1) 

in which f = r^cr'^, yU = 1,2. 

We want to use Eq. (|3.4|) with n = 1,2, that we rewrite here (neglecting 
the trivially zero contributions and changing slightly notations): 

a^Ci(r) = hm { < [: a„(xi) : (V^(r)V^(O) - C"^-0™)] >} (5.2) 



/Ri rR2 — 

d^x, / d^x2 < [: ouxi) ■■■■ 0UX2) ■■ mr)m 
J 

-C - C^'^^)] > - f'd^Xid^C''- < [: Om{x,) :: 0^ :] > -{R2 ^ Ri) 

(5.3) 

where Om and $2 = (^r^?/')(?/'r'^?/^) are the marginal scalar operators of the 
theory [V^ = (1, 7^, Cs)). Notice that lower indices of Wilson coefficients have 
not been reported. 

We consider this example very instructive because the values of the (dif- 
ferentiated) zero mass Wilson coefficients C"^'",C*2,c}mC"^™ will be fixed im- 
posing cancellation of I.R. divergences (of course the first two can be obtained 
in independent way from the massless theory). 

In principle we should give proper normalization conditions for the one 
and two insertions of the operator Om = ip'ip in such a way that the Action 
Principle (Eq. (|4.6| ) with k = 1 and $a = Om) holds. The point is that 
Action Principle is automatically satisfied in the massive theory if the usual 
MS definition for composite operators is adopted (see e.g. |^). If we use for 
a moment the knowledge of the massive theory propagator, we obtain after 
standard calculations: 

d'^p 
(2^ 

1 /"i / + x{l — x)p'^ , 



< [: Om{r) :: 0^(0) :] >„= j ^e~^P^F{p,m) (5.4) 



F{p, m) = -— [ (ia; ( 7ij + 3 + log 
2n Jo \ 



2n Jo y^' ' Ann'^ J 

2 

< [Om] >m= -^hE + 1 + l0g(^)) (5.5) 

It is easy to check that 

dm < [Om] >rr.= jdh < [. Om{r) :: Om{Q) :] >m= F(0,m), (5.6) 
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i.e. Action Principle holds in this case without addition of contact terms. 

Having we used only the knowledge of the massless theory, we would get 
Eq.fpD with 

1 

— (7i^+l + log(-^^ 



F{p) ^ F{p, 0) = - — (7^ + 1 + log(^)). (5.7) 



We could also get a partial information on < [Om] >m imposing that the fj, 
dependence (i.e. the U.V. behavior or in other words the anomalous dimen- 
sions) is the same in massive and massless theory: 

fxd^dm < [O^m] >m= f^O^ fd?x < [. 0^{x) :: 0^(0) :] > (5.8) 

(notice the m dependence of left side). The cutoff integral can be computed 
exactly by using the properties of Bessel functions : 



fd'x < [: Om{x) :: O„(0) :] >= J0-^F{p)^MpR) 



= i-(log(7r/i2i?2)+7^-l)) (5.9) 

Integrating the differential equation ( |5.8|) and imposing the dependence on 

mj [i we get 

1 Tfl 

dm < [Cm] >m= log h COUSt, (5.10) 

vr /i 

from which, integrating again, we get 

1 Tfl 

< [Om] >m= m(log - + Co J. (5.11) 

vr /i 

The constant Cq^ can be fixed only imposing the Action Principle on the 
deformed theory (this gives the correct V.E.V., Eq. ( |5.5|) ) and its existence 
has been anticipated in Section ^ (notice that the two point correlator has 
zero I.R. degree of divergence and we expected I.R. problems in this case). 

We computed some of the required integrals passing to complex coordi- 
nates r = r^+ir"^ ■ ■ ■ and using Stokes' Theorem to reduce surface integrals to 



line integrals. In [Appendix A| we give an explicit example of such calculations 



The results (up to not relevant subleadings terms in the limit R2 » 
Ri » \r\) are: 

""Vxi < [: OUxi) : ^(r)^(O)] > 7^ log ^ 

ZTT r 
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/ < [: OUxi) ■■■■ Omix2) ■■ (^^(r)V(O) - C\r))] > 



112 112 
?2 



For what concerns the remaining integrals, it suffices to note that 

/Ri rR2 
J <[:0^: (xi) : : (x2)$^(0)] >^ log fiR, log fiR2 + ■■ ■ (5.13) 

as it is easy to check. 

Equating terms of same behavior in R, by Eq. ( |5.9D , Eq. ( |5.13| ) , Eq. ( [5.12[ ) , 
we get from Eq. 



C^- = -1/2 d^C' = ^ (log(7r/i>n +7^-1) (5.14) 



An 

and from Eq. (|5.3|) 

a^c^- = If C*2 = 

die = (log7r/i>|2 - 2 + 7s) f. (5.15) 

Combining these results with the VEV of Om Eq.( |5.11|) (and with the off 
critical information ([5.5D), we get 



1 771 V p'^^ 1 

<ij{r)m>m - --— (l + 2"^'kP(log(^^)-- 

lE 

,3 



1 


f 


2tx 




m 


log( 




+2^ 





^ )+0(m^) (5.16) 

that is the correct result as can be checked from the properties of K functions. 
Another check can be obtained from the constraints 

^^j-^^lC^ + nar'C^'"/i|^ < [Om] >m= 

^il^^^""^ = (5.17) 
obtained from the free theory relation 

° " ^ "^"^ ^'"^ ^^^^'^""^ ^ C7^'"(m) < Om >m) (5.18) 
equating terms of same order in m and logm. 
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5.2 Application to Ising Model 



We want to reconstruct here the mass corrections to short distance behavior 
of the two points spin-spin correlation function of the Ising Model by using 
the previously introduced method to compute Wilson coefficients. Similar 



results have been already obtained in a previous article |T0[; we only want 
to treat this model as a concrete example of our reformulation to be com- 
pared with the preexisting approach. To put in evidence the renormalization 
scheme independence of our treatment (provided Hypotheses 1-3 are satis- 
fied), we will use here the Hadamard finite part (see [Appendix Bp to deal 
with U.V. divergence and the I.R. regulator Q{x) = Jo(|Q||x|) (see Eg. (|3.12|) ) 
that is equivalent to insert G(x) = e**^^ and then average over directions of 
before taking the limit Q 0. 

It is a well known fact that the Ising model in proximity of its critical 
point is described by a theory of free Majorana fermions with mass m and 
euclidean action 

1 ^2. 



S= — jdx'iljd4) + Tpd'ilj + imi)ip. (5.19) 

ZTT 



Moreover the exact correlator [R3] is known from the scaling limit of the 



lattice theory and can be used check. 

The conformal theory at fixed point m = is described by the primary 
operators l,cr, £^ =: -^ipip '■ of dimension x = 0,1/8,1 from which the sec- 
ondary operators can be built up by repeated application of the Virasoro 
generators: 

^{n}{n} ^ . . . L_^^L_n, ■ ■ ■ L-n,s^, (5.20) 

where $ is any of the primary operators and the dimension of the secondary 
operator is a; = + '^j + S ^j- Away from the fixed point the expression 
of the operator conjugate to will be of the form [Om] = < [Om] > 1- 

At first order we have 

dr^Clir) = Hm < [-^(x)(a(r)a(0) - Cf,(r)f (0))] > (5.21) 

(the average over directions of to be taken after integrations is not re- 
ported for simplicity of notations). 

Before proceeding to the computation of the integrals a precise definition 
of the operator Om should be given in order to obtain well defined distri- 
butions from its insertions, as discussed in Section ^. We define the double 
insertion of Om by 

2 

< [Om{r)Omm >m= Fp^— (iT^^ (mr) - K^imr)) (5.22) 
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where Fp^ means the distribution obtained taking the Hadamard finite part 
of the U.V. singularity. (We used some additional knowledge of the massive 
theory: the massless correlator can be derived by taking the limit m — > of 
previous expression). 

Proceeding as in Section ^ and imposing the Action Principle 

/I jYI 
dh < [Om{r)0^m >m= - — (l0g( — ) + 1-^(1)) 

(5.23) 

(see ([B.6|) for computation of the integral), we derive the V.E.V.: 



m , ,me 



IE 



< [O^] >ra= -;^log(^^). (5.24) 



We recall that (as explained in Section ^ and in Section |5.1| ) the numerical 
constant inside the logarithmic term cannot be fixed only by knowledge of 
the m = theory. 

We can now come back to the computation of integrals in Eq. (|5.21| ). 

The (locally integrable) correlation functions involved in the first integral 
is exactly known from the conformal theory |]33| : 

3/4 

< ^(ri)(T(r)a(0) >= — — (5.25) 

47r|ri — r\\r\\ 

and the integral can be computed (in the limit Q ^ ) by standard Feynman 
diagram techniques and properties of Bessel Functions [^ : 



\x-r\\r\ Jo (t(l-t))V2 7o s 

= 2 J^dt ^^^^_^^^^^^ Ko{^t{l-t)\Q\\r\) ^ -27rlog ''^"g' (5.26) 

For what concerns the second integral in ( ^.21| ) we have (after the limit 
m — s> is taken) 

Jd^xe'^'' < [S{x)S{0)] > = y^W^^Fp^^-^ 

= ^(logf + 7.) (5.27) 



see 



3^ for the Fourier transform of Hadamard finite part). Substituting 
dlD-O) in Eq. (^:2TD we get 

dn.Cl = K'/\\og ^j^ + jE- log 4) - Cl^ilog M + 7i.). (5.28) 
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Cancellation of the log Q factor fixes 

C5^(r) = Trr^/^ (5.29) 
while equating the regular term we get 

d^Cl = r'"\\og{^) (5.30) 



From ( |5.30D we can thus reconstruct at first order in m the deviation from 



the scaling behavior of the fixed point theory for the Wilson coefficient 

ClAr, m;f^) = ^,{l + ^ log(^) + 0((mr)2)) (5.31) 



and (using also ( |5.29|) ) of the short distance behavior of the two point func- 
tion: 

< a{r)(x{0) ~ Cl^{r, m; ^) + mdmCZ{r, m = 0) < 0„ >^ 

-(1 + — log(— e^-) + 0{{mrf)). (5.32) 



rV4 ^ 2 



These results are in agreement with those of [|10l and with the exact ones of 



34| 



6 Conclusions 

We presented a consistent method to deal with I.R. singularities of expan- 
sions in a relevant parameter around a field theory without super-renormalized 
couplings. Our approach is a reformulation and an extension of other recent 
works p, [l^]- The goal is the reconstruction of the short distance behavior 



(Wilson coefficients) of the deformed theory obtained by adding a relevant 
operator to the action. The use of Action Principle and (weak) convergence 
property of O.P.E. give closed and I.R. finite expressions for the derivatives 
with respect to generalized masses of Wilson coefficients that, can thus be 
reconstructed by means of Taylor expansion (that is convergent if analyticity 
of Wilson coefficients can be assumed, otherwise asymptotic). 

Our contribution is an orthodox treatment of U.V. divergences of com- 
posite operators that, as a consequence, allows for an explicit and simple 
inductive proof of I.R. finitess and of general equation Eq. (|3.12|) . We put in 
evidence that Eq.( |3.12D does not refer to a particular U.V. renormalization 
scheme (provided Hyp. 1-3 are satisfied) and holds for general I.R. regulators. 
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We discussed also the presence of well defined but unknown constants 
related to validity of renormalized Action Principle that parameterize the 
long distance behavior of deformed theory (essentially operators' V.E.V.) 
and cannot be computed by using the knowledge of the massless theory 
only. Nevertheless these constant are universally present in short distance 
expansions of different products of composite operators and are always in 
finite number at any finite order in the masses. The possibility of fixing 
these parameters (e.g. by some symmetry relation between operators in the 
deformed theory or the knowledge of some particular deformed correlator, as 
in the Ising model of Section p.2[ ) depends on particular models and should 
be discussed case by case. 

In spite of this problem, we think that its simplicity and generality, com- 
pared to many specific approaches, make of this method a powerful and 
promising tool. 

We remind in particular that the (three) main hypotheses on which it 
relies have been found to be satisfied for the bulk of the interesting conformal 
field theories if adequate renormalization schemes are used: it thus seem that 
the technique could be helpful in the investigation of relevant deformations 
of those theories. We hope to come back to the problem soon. 

Acknowledgements: The authors are grateful to CM. Becchi for sug- 
gesting to them the problem and for many useful discussions. We acknowl- 
edge also for interesting discussions R. Collina, G. Delfino, K. Konishi, C. 
Imbimbo, Hiroshi Suzuki, J.-B. Zuber. Part of the work of one author (R.G.) 
was realized during visits to S.Ph.T. of Saclay and C.E.R.N., whose staffs 
are acknowledged for the warm hospitality. 



Appendix A Complex integrals 

We give here an example of the computations required to obtain Eqs. ( |5.23| ) 
the detailed calculation of the integral 

d^z- ^ (A.l) 

\z\<R [z — \r\)* z — z' 

in which 2;', r are complex variables (2;' 7^ l^^l) and (P'z = ^dz A dz*. A trick 
to simplify things is to use Stokes' theorem: given a one form Q regular in 
the domain S one has 

[dn= I n. (A.2) 
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If we pose Q = uj{z, z*)dz + uj{z, z*)dz* we can fix VL in our case requiring: 

,du) duo , i 1 1 



(A.3) 



dz dz* 2 {z — \r\)* z — z' 
With the choice a; = and = — | \og{z* — \r\) / {z — z') we get 

/ = Jhgiz* - \r\)/{z - z')dz, (A.4) 

in which the path 7 should not enclose singularities (pole in and cut from 
|r| that we choose along the real positive axis) to allow for application of 
Stokes' Theorem. We split 7 = 71 U 72 U 73 U 74 where (in terms of z): 



7i 


^{z/z 


= z' + ee' 


^0G (27r,0)} 


72 


^{z/z 


= pe [\r\ 


,R] or ^ = e2>,pe(i?,|r|)} 


73 


^{z/z 


= Re'''', 


e (0,27r)} 



(A.5) 

(notice the orientations of the intervals). The contribution of the path 74 
turning around \r\ and connecting the two branches of 72 is easily seen to be 
zero. 

Performing easy integrations we get the desired result: 

I \ 7^ — 1 7^ I I ^ \ I T I \ 7^ \ 



All the integrals in Eqs. (|5.23|) can be computed at relevant order by (re- 



peated) application of Stokes' Theorem. 



Appendix B Computation of the Hadamard finite part 



In one dimension Hadamard finite part of a non locally summable 
(in X = 0) function 



X 



x] 



X 



(B.l) 



(in which Xi, is a finite set of complex indices with ReAj, > 1 but different 
from 1 and h{x) is regular) is 



Fp J^g{x)dx = \imljj{x) - /(e) 



(B.2) 
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in which the infinite part is defined as 



/(e)-ET^(V-^ + ^ilog-- (B.3) 

„ A,, — i e e 



In this case we obtain 



Fp [\ix)dx = -y T^^fi)^""^ - Ai log \ + f\(x)dx. (B.4) 
JO ^ X,^ — 1 b b Jo 

One can extend a locally non summable function as a tempered distri- 
bution by taking the Hadamard finite part of the integral of the function 
times the test function. For our purposes we are involved only in radially 
symmetric functions and this case is a trivial generalization of one dimen- 
sional situation. Extension to higher number of dimensions can be found in 
pSi and [Q. In particular in is given a detailed treatement of distri- 
butional extensions of homogeneous functions of many variables that could 
be useful in the treatment of conformal field theories singularities. Notice 
also that Hadamard finite part, subtracting only divergent part is a Minimal 
Subtraction in the sense explained in Section ||. 

We compute now the integral 



47r2 



2 

{K^{mr) - K^{mr)) (B.5) 



The integral of this distribution is easily done by taking the finite part at /ir = 
i.e. by subtracting the divergent part in /xr = (we added a renormalization 
point scale fi to keep correct dimensions): 

Fp, / d\—[Kl[mT) - Klimr)) = — Fp, / dx x{KI{x) - KI{x)) 



2 



oo 



/e<|r| 

= ^Fp,(^(2ir^(x) - K,{xf - K,{x)K,{x)))^^ 

1 Tflf 

= -:^Fp,(log(— ) + l-V^(l)) 

= -^(log(^) + 1-^(1)) (B.6) 

where the properties of the Ky{z) functions (in particular their behavior in 
z = Q and m. z = oo) have been used, see for instance 
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